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ABSTRACT 



A digital halftoning method and apparatus with error 
diffusion reduces visibly discernible artifacts. In the 
error diffusion method and apparatus, a range function 
provides random factors which arc used to spread crros 
to neighboring pixels. 

6 Claims, 6 Drawing Sheets 
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instance, can be assigned to pixel locations where de- 
PROCESS FOR PROVIDING DIGITAL HALFTONE tected gray-scale values are less than 1 27, and the binary 
IMAGES WITH RANDOM ERROR DIFFUSION number "0" can be assigned to pixel locations where 

detected gray-scale values are 127 or more. 
This application is a continuation-in-part of co-pend- 5 As an example of binary coding for error diffusion in 
ing U.S. patent application Ser. No. 07/257,843 filed digital halftoning, a pUe! location with a detected gray- 
Oct. 14, 1988 and U.S. patent application Ser. No. scale value of 100 would be coded as a binary "1**. 
07/421,931 filed Oct. 16, 1989. Similarly, a pixel location with a gray-scale value of 145 

would be coded as a binary "0". In monochrome print- 

BACKGROUND OF THE INVENTION jq j^g systems, a pixel location encoded as a binary "1" 
The present invention generally relates to a method generally is blact and a pixel encoded as a binary "O" 
for producing images by digital halftoning and, more generally is white. 

particularly, to a production method based upon digital The above-described coding techniques for digital 
halftoning with error diffusion. halftoning with error diffusion inherently result in cod- 

13 ing errors for each pixel location having a gray-scale 
STATE OF THE ART value other than zero or 255. For instance, the magni- 

In conventional halftone printing, picture elements, tude of the coding error resulting from binary encoding 
or pixels, of a continuous-tone image are converted to a pixel location having a gray-scale value of 100 would 
dots whose diameters vary in size according to the be 100 units as measured by gray-scale values. Similarly, 
amplitude of light pcnetratmg a screen. More particu- 20 the magnitude of the coding error for binary encoding 
larly, the dots are larger where more light penetrates a of a pixel location having a gray-scale value of 110 
screen and are smaller where less light is present. Thus, would be 1 10 units as measured by gray-scale values, 
conventional halftone processes can be described as The basic concept underlying digital halftone print- 
being of the analog type since in the processes, dot size ing with error diffusion is to distribute, or diffuse, cod- 
is the analog of color intensity for individual picture 25 ing errors from the locations at which the errors arise, 
elements of an original image. This concept is further discussed in an article, entitled 

To produce multi-color images from continuous-tone "An Adaptive Algorithm for Spatial Greyscale", R. W. 
originals by conventional halftone printing, a color-fil- Floyd and L. Steinberg. Proc, SID, vol. 17/2, pp. 75-77 
tcred screen is prepared for each of four primary colors. (1976). Various algorithms have been proposed to ac- 
Then the colors are printed sequentially while maintain- 30 complish error diffusion. 

ing precise mechanical registration of the screens. Al- Although images produced by digital halftone print- 
though conventional halftone printing can produce ing with conventional error diffusion techniques may 
excellent results, the processes involve economies-to- not have a grainy appearance, the printed images often 
scale and are relatively expensive for short production contain discernible patterns. Such patterns are often 

35 referred to as "artifacts." Artifacts sometimes appear as 

Digital halftone printing, in contrast to conventional straight lines and at other times as wiggly or worm-like 
halftone printing, is better suited for shorter production shapes. At still other times, artifacU have a random, 
runs. In digital halftone printing, the initial step is to pebble-like appearance. Regardless of their particular 
scan an image with an optical scanner that detects light form or shape, however, artifacts detract from the ap- 
intensity values for selected colors. Then for each pixel 40 pcarance of printed images and therefore limit the com- 
location. the detected data are converted to gray-scale merdal usefulness of digital halftoning, 
values to provide a basis for image, production. Typi- SUMMARY OF THE INVENTION 

cally, gray-scale values range from zero to 255. 

More particularly, in digital halftone printing, dots In accordance with the preferred embodiments, gen- 
arc printed in areas corresponding in scale to the origi- 45 crally speaking, the present invention provides an im- 
nal pixel locations. The printed dots normaUy arc of proved error diffusion process for ehminating visually- 
fixed size, but their placement and density within the perceptible artifacts in images produced by digital half- 
printed areas are vimed in accordance with the detected tone imaging. These techniques are as applicable to 
gray-scale values. Thus, in images that have Keen pro- printing as they are to desk-top publishing, and to any 
duced by digital halftone printing, dot density within 50 other application that uses a pixel pattern of an image, 
each printed area determines the spatial resolution and regardless of the origin of the pixel pattern or the me- 
perccived grayness of the area to an observer. dhun where the image is displayed. 

One shortcoming of conventional digital halftone In the first preferred embodiment, the present inven- 
printing is that printed images often have a grainy or tion provides a process with error diffusion, comprising 
"noisy" appearance. To reduce the grainy appearance 55 the steps of: 

of digital halftone images and to improve detail resolu- detecting gray-scale values at preselected pUel loca- 
tion, it has been suggested to use error diflftision tech- tions in an image; 

niques in conjunction with digital halftone printing. For for each detected gray-scale value at the preselected 
example, the text Digitai Halftoning, by Robert Ulich- pixel locations, determining upper and lower limits of a 
ney (MIT Press, 1987), pp. 265-268. describes error 60 range fimction which depends on the detected gray- 
diffusion techniques with perturbation. Those particular scale values, and selecting first and second values which 
error diffusion techniques are sometimes referred to as are randomly located between the upper and lower 
spatial dithering. limits of the range function; , ^ . , 

To apply error diffusion techniques to digital halftone if the detected gray-scale value at a preselected pud 
printing, grayscale values at detected pixel locations 65 location is less than a predetermined threshold value, 
are binary coded. Typically, the binary coding thresh- encoding the pixel location as a binary 1, assigninq a 
old is 127. i,e.. halfway between the minimum and maxi- first error value for the preselected pixel location equal 
mum gray-scale values. The binary number "P, for to its detected gray-scale value multiplied by the first of 
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the selected randomly located values, assigning a sec- for selected pixel neighbors of the ith pixel in an 
end error value for the preselected pixel location equal image, diffusing the calculated error value, E/, of the ith 
to its detected gray-scale value multiplied by the second pixel location to its selected neighbors as follows: 
of the selected randomly located values, and assigning a 
third error value for the preselected pixel location equal 5 ey=(Hy)H/ 

to its detected gray-scale value multiplied by 1 minus ^.».., .v 

the sum of the first of the selected randomly located where e; is the error value diffused to the jth pixel neigh- 
values and the second of the selected randomly located bor from the ith pixel location, where the jth puel 
y^ygj. neighbor is randomly determined from a predetenmned 

if the detected gray-scalc value at a preselected pixel 10 set of pixels neighboring the ith pixel, and where w;is a 
location is greater than the threshold value, encoding weighting factor determined for the jth pixel location; 
the pixel location as a binary 0, assigning a first error continuing the error diffusion process for the other 
value for the preselected pixel location which equals the pixel neighbors of the ith pixel location; 
first of the selected randomly located values multiplied continuing the error diffusion process for the Q + l)th 
by a quantity equal to the detected gray-scale value less ^5 pixel location and so forth, such that coding errors are 
255, assigning a second error value for the preselected diffused from substantially all of the pixel locations m a 
pixel location which equals the second of the selected scanned image; and . . j 

randomly located values multiplied by a quantity equal printing images by digital halftone pnntmg based 
to the detected gray-scale value less 255, and assigning upon intensity values of pixel locations that have been 
t third error value for the preselected pixel location 20 modified by error diffusion. 

equal to its detected gray-scale value less 255 multiplied Generally speaking, it is preferred that the predeter- 
by I minus the sum ofthe first of the selected randomly mined set of neighbor jpixels comprises pixels that are 
located values and the second of the selected randomly immediately adjacent neighbors to a selected ith pixel 
located values; location. More particularly, it is preferred that the pre- 

then, for each preselected pixel location, distributing 25 determined set of neighbor pixels comprises the pixel 
the first, second and third assigned error values to three that immediately follows a selected ith pixel location, a 
adjacent preselected pixel locations, one of which lies pixel inunediately below the selected ith pixel location, 
on the same line as the preselected pixel location and the and the pixels that immediately precede and follow the 
other two of which are on an adjacent line; and pixel immediately below the selected ith pixel location. 

printing images based upon values of the preselected Further in the preferred embodiment, the weights for 
pixel locations that have been encoded after distribution the jth and subsequent pixel locations are determined as 
of the assigned error values. follows: 

More particularly in the first preferred embodiment, a) for the first randomly chosen adjacent neighbor 
the lower limit of the range function is between about pixel, assigning a weight w<// which is randomly selected 
zero and 0.2 and varies linearly with the detected gray- from a distribution of numbers between 0.0 and 1.0; 
scale value of a preselected pixel location. More partic- b) for the second chosen adjacent neighbor pixel, 
ularly, the lower limit, Wmin> of the range function pref- assigning a weight w<^ which is randomly selected from 
erably is determined by the following equation: a distribution of numbers between 0.0 and (1.0— wji); 

^ c) for the third randomly chosen adjacent neighbor 
|2 _ n pixel, assigning a weight Wrfjt which is randomly se- 

H'mifl = — f — X 0.2 lected from a distribution of numbers between 0.0 and 

(\.0-v/di-wdj); and 

and the upper limit, Wmax, of the range function is deter- d) for the last remaining pixel of said set, assigning a 
mined by the following equation: ^5 weight w^/ to it where 

I f — xa6 

BRIEF DESCRIPTION OF THE DRAWING . 
where z is the detected gray-scale value at a preselected 50 ^ 

p«dlocation.Alten,ativdy.thettpperandlowerlinut8 The present mv«ttton can be further undm^ 
Sfthe range ftmction can be deteimined from the graph f ^^ce to &e foBowong descnption «idatt.ched 

of FIG 3 drawmg which dlustratcs the preferred embodiments. 

In the second preferred embodiment, the error diffu- In Redrawing. r . ,r 

son process includes the steps of: 55 FIG- ^ ^ » fiinctional block diagram of a digital half- 

scuming an image to detect intensity values at se- system which operates accordmg to the process 

lected pixel locations- ^ present invention; 

encoding selected ones of the pixel locations as one FIG- 2 schematically shows a field of pixels and is 
binary value when the detected intensity value at a offered as an aid to the description of the enor diffusion 
location is less than a predetermined threshold value, 60 process ofthe present invention; 
and calculating error values for such pixel locations that FIG. 3 shows a function for use m the process of the 
equal the detected intensity values; Present invention; ^ ^ , ^. 

encoding selected other ones of the pixel locations as FIG. 4 is a more detailed block diagram ofthe system 
a second bmary value when the detected intensity val- of FIG. 1; • , . 

ues at those locations are greater than the predeter- 65 FIG. 5 is a discrete block diagram implementation of 
mined threshold value, and calculating error values for the first embodiment of the present invention; 
those locations that equal the detected intensity values FIG. 6 is a discrete block diagram unplementation of 
less the upper limit of the intensity scale; 'he second embodiment of the present mvention; and 
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FIG. 7 is a block diagram of the blocking circuit for the pixel designated by arrow da can be described as 
FIG. 6 to account for edge conditions. having coordinates (x,y+ 1), and, finally, the pixel des- 

ignated by arrow d4 can be described as having coordi- 
DESCRIPTION OF THE PREFERRED ^^^^ (x- l,y-H 1). The set of directions d, through da 

EMBODIMENTS 5 through d4 are herein called the "direction sets" 

FIG. 1 generally shows a system for digital halftone with di through da for the first embodiment and d\ 
printing. The components of the system comprise an through d4 for the second embodiment, 
optical image scanner 5, a general purpose digital com- FIG, 4 is an expanded block diagram over the block 
puter 7, and a digital printer 9* The illustrated compo- diagram shown in FIG. 1. FIG. 4 shows several addi- 
nents are conventional* but by employing the method of 10 tional input/output devices which might be used for the 
the present invention, the components operate in the various possible applications of the present invention 
system to eliminate the appearance of artifacts in images including image production, printing, desk-top publish- 
printed by digital halftoning with error diffusion. ing and CRT (cathode ray tube) display, to name just a 

InnormalopcrationofthesystemofFlG.l, scanner few. In addition to optical scanner 5 and printer 9 
5 makes color separations of an image from a medium 13 IS (which are both shown here and in FIG. 1), there is a 
such as a photographic transparency, a photographic magnetic or optical disk drive 28, a keyboard 30, a 
print, or a high-resolution printed original. For each of monitor 32 and a tape drive 31. Also, FIG. 4 shows 
the separated colors, the picture elements (i.e., pixels) of internal details of digital computer 7, namely, a proces- 
the original image are digitized to gray-scale values. sor 20 which communicates via a data bus with RAM 
Then the digitized values are provided to digital com- 20 26, ROM 24 and input/output interface 22. Input/out- 
puter 7. Within computer 7, the digital information can put interface 22, in turn, then might communicate with 
be processed with conventional graphics software so disk drive 28, keyboard 30, scanner 5, monitor 32, tape 
that images can be retouched, adjusted, recomposed, drive 31, and printer 9 via a secondary data bus that is 
and so forth. Then the processed images are provided to external to computer 7. 

digital printer 9 for printing onto a sheet 15 by the tech- 25 3y including the additional input/output devices the 
niques of conventional digital halftone printing. ^ta, or image, that k to be enhanced by the present 

When the system of FIG. 1 is operated in accordance invention can be other than a photographic image 
with the present invention, the digital image informa- which is input into the system by first scanning a photo- 
tion that is provided to computer 7 is also binary en- graph, as discussed above. For example, the image to be 
coded and then is operated upon according to the half- 30 pri)cessed may be a computer-generated image that is 
toning process of the present invention to provide error ^j^g displayed on monitor 32 where the processed 
diffusion prior to printing. The object of the error diffu- image may cither be displayed on monitor 32, printed 
sion process of the present invention is to create rcpro- printer 9, or it may be represented by a previously 

duciblc images that provide the illusion of a continuous- stored pixel pattern which is stored magnetically or 
tone image without discernible artifacts. 35 optically on a disk. There are many options which are 

An error diffusion process according to the present discussed further elsewhere herein, 
invention will now be described with reference to FIG. purposes of discussion, the error diffusion pro- 

2. In FIG. 2, each pixel 11 in the field should be under- of both embodiments of the present invention will 

stood to represent one pixel location on an image which described as beginning at pixel P in FIG. 2. Further 
is presented for optical scanning. For purposes of dis- 40 purposes of discussion, the magnitude of the de- 
cussion, it can be assumed that scanning proceeds from tectcd gray-scale value of pixel P will be described as an 
left to right across each pixel row and sequentially, row arbitrary number * V* which is between zero and 255, 
by row, from the top to the bottom of the field. Such a inclusive. Still further, the coding error for any pixel 
scanning pattern can be described as a raster-type scan- location in the field will be referred to as Exj> where the 
ning pattern. Other scanning patterns, however, might 45 subscripts identify the cartesian coordinates of the pixel, 
be used; for example, scanning can proceed from left to When a pixel at a given location (x,y) has a value z 
right across one row of pixels, then in the opposite that lies between a threshold value T and the maximum 
direction across the next lower row, and so forth in a gray-scale value (e.g., 255), the magnitude of the encod- 
serpentine pattern. jog error Exjf will be a negative number equal to z minus 

The location of each pixel 11 in the field in HG. 2 can 50 255 gray-scale units; that is, 
be described by cartesian coordinates; accordingly, an 

arbitrary pixel '*P*' is designated as haying location Sg^^t-lSS tot z^T, 

(x,y). Because pixel P is in the interior of the field, it has 

eight immediately-adjacent pixel neighbors. When the When the pixel value z is a number between zero and 
pixel field is scanned in a raster*type manner, four of the 55 the threshold value, the magnitude of the encoding 
immediately-adjacent neighboring pixels will be error is simply equal to the value z in gray-scale units; 
scanned before pnel P, and four will be scanned after that is, 
pixel P. 

In FIG. 2, the pixel to the immediate right of pixel P Ex^^x for x<T. 

is designated by arrow dl, the pixel diagonally below 60 . ^ , ^ . « 

pixel P in the scanning direction is designated by arrow The selected threshold value T is usually a number 
d2, the pixel directly below pixel P is designated by about midway along the selected gray-scale and typi- 
arrow dj, and the pUel diagonally below pixel P in the cally cquaU 128. Although the threshold value usuaUy 
direction opposite the scannmg direction is designated is held constant, it could be randomly vaned, or dith- 
by the arrow d4. TTius, the pixel designated by arrow di 65 crcd, in a relatively narrow range about a chosen cen- 
can be described as having cartesian coordinates tral value (e.g.. 128). In that case, further randomness 

' would be introduced into the error diffusion process 



(x+ l,y), and the pixel designated by arrow d2 can be 
described as having coordinates (x-|- l,y-l- 1). Likewise, 



that will be described below. 
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- . scale value for that pixel location plus the propagated 

First Embodiment ^^^^^ Likewise, the gray-scale value of the pUel 

In the first embodiment of the present invention, the located at (x+ l,y H- 1) becomes the detected gray-scale 
gray-scale values z of pixels are used as the arguments value for that pixel plus the propagated error erri. Fi- 
of range functions. One example of range functions, 5 nally, the gray-scale value of the pixel located at 
Wm/rt and Wffjojri is shown in FIG. 3. In this example, the (x,y-i- 1) becomes the detected gray-scale value at that 
range functions provide variable minimum and maxi- pixel plus the propagated error errs, 
mum limits which are linearly dependent upon the de- It should be understood that the above-described 
tected gray-scale values z of pixels. For pixel values error difhision process continues with successively 
ranging between zero and 128 , the smallest minim u m lo scanned pixels. Thus, the binary encoding error for the 
value ofthe range function Wm/nis zero, and the largest pixel located at (x+l,y) will reflect the propagated 
minimum value is 0.2. Also, for pixel values z between ^rrcr crri as well as the crron propagated from other 
128 and 255, the smallest minimum value of the range previously scanned pixel locations. In turn, the error 
function w„/n is zero, and the largest minimum value is value for pixel location (%+ l,y) is weighted by the 
0.2. Algebraically, the range function Wmfa can be ap- 13 pseudo-random process and then diffused to the pixels 
proximated as foUows: located at (x-f 2,y), (x-f l,y-l- 1). and (x4-2,y + 1). The 

process continues for each pixel in a row and then be- 
\z~T\ _ gins again at the first pixel on the next lower row. 

wm/fl ^ J X In practice, special error assignment rules must be 

20 provided for the first and the last pixel locations in a 
Further in FIG. 3, for pixel values zbetween zero and row. For the last pixel location scanned in a row, a 
128 , the smallest maximum value of the range function convenient rule is that an encoding error is assigned 
Wmax is 0.4, and the largest maximum value is 1.0. For only to the pixel located immediately below. For the 
pixel values between 128 and 255, the smallest maxi- grst pixel scanned in a row, a convenient rule is that its 
mum value of the range function Wmax is 0.4, and the 25 assigned error value reflects only the error assigned 
largest maximum value is 1.0. Algebraically, the range |^om the pixel immediately above, 
function Wmax can be approximated as follows: The above-described error diffusion process can be 

summarized as follows. First, in accordance with the 
\z~ T\ range distribution function of FIG. 3, the range of 

^max'^i f — X • 30 wm//i,wmai is determined by the gray-scale value of the 

pixel located at (x,y). Then, the following calculations 
According to the present invention, the range f\mc- ^rc made: 
tions v/min and Wm« are employed to calculate pseudo- 1^ n =RANDOM [wm;fl(pixcl(x,y)), w„ax(piJtel(xj,))] 
random weights ri, rj and rj. For calculation of pseudo- 2. ra^RANDOM [wm/n(pixcl(x,y)), Wmax(pixel(xj-))i 
random weight n, a random number is selected, then 35 3 rz^hO—r\~n 

the random number is multiplied by the numerical range 4 p^^ e^ch input pixel location (x.y), if the detected 
between Wmfn and Wmax. and the resultant product is gray-scale value of the pixel at that location is less 
added to Wmin- In practice, the random selection pro- preselected threshold value, then 

cess does not employ a random number generator but, pixcl(xj») is black, and 
to minimize processing time, uses a look-up Uble Of 40 error=value of pUelz-y) 
uniformly distributed random numbers expressed as 

decimals. An identical procedure is followed to calcu- pixela ) is white, and 
late weight rj. Then the sum of the pseudo-random crror^pixel(jtj)-255. 

weights ri and rs is subtracted from the unity to yield a 5 ^j^^ ^jistnbmed* error values are calculated as 
third pseudo-random weight n; that is, 17= 1.0— n—n- 45 foUows: 

The pseudo-random weights ri, r: and ra are used to erri ^ri * u, 
propagate, or spread, error values to the pixels at loca- em=ri • ExZ and 
tions (x+l.y). (x+l,y+l) and (x,y+l). respectively. erra^rs* Ejci 

More particularly the error propagated to the pacl ^ ^^^^ distributed to neighboring 

located at (x-f l,y) from the pixel located at {^y)^ 50 ^^^^^ ^ 

determmed by multiplymg the weight n by the encod- pUeWi rt=pixeU+i^)+em 
ing error for pixel P; that is. em=ri • E;^ Likewise, *^Ln-^Sd^ll,^^^^ 
the error propagated to the pixel located at (x.y +.1) is 

m2=r2»E:^.Finally,theerrorpro^^^ ^d^Tsho.^ a b£^k diagLn implementation of 

locat«Jat(x+l.y-f l)iserr3=r3»P,o.Tlius. Aep^^ 35 ^^ ^^^ 7 fo, ^ embodiment error difftision 
agat«j errors erri. em, errj can be understood to be t^^niq^e of the present invention. This implementation 
pseudo-randomized error values which are diffused or ^ for a pixel pattern that is scanned, displayed or stored, 
distributed to pixels neighbonng pixel P m ^^^^^on ^ ^ j^,Q ^ ^ j;^^^ 

of the arrows d,, dj, and dj, respectively, of FIG. 2. ^ ^ ^ ^ ^ ^ ^^^^ ^ 

It should be noted that the i,ropaptede«^^ 60 g^y.scale values for each pixel are sequentially applied 

either positive or negauve. Also, it should be under- ^ ^ ^4- 1 byte-long register 32 from I/O interface 22 
stood that pixels with, detected gray-scrfe values of ^ ^^^^ fo^ p 

citherz^oor 255 do n^generateerror values. Thatis, p with the z value for P<x+i^+i) bcSg the next 
error vducs are npt diffu^^ from pixels havmg gray- ^ ^^^^^ received by register ^TFrom register 32 the 
scide values ofzcro and 255. but error values can be 65 ^ values of each of the four pixels needed for the opera- 
diffused to such pixels. 1 , , tionofthe first embodiment ofthe present invention arc 

For bmary encodmg purposes, the gray-scale value of available 
the pixel located at (x -H 1 ,y) becomes the detected gray- 
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There are two paths in this operation which come signal, and when x equals N or y equals M, to block the 
together at the end for the final calculations. First, there changes to the pixel values in register 32 from delays 96 
is the path for calculating the error associated with the and 98 (x=N) and delays 94 and 96 (y=M), respec- 
operational pixel, Px^. Also. Zx^is applied to compara- tively. The blocking discussed above is necessary to 
tor 46, the positive terminal of subtracter 48» and multi- 5 prevent the value of the pixels in the last column from 
plcxer 50. In sublractor 4S the value 255 is subtracted effecting the pixels in the first column (x=N). and from 
from Zxj^ and the resulting value from subtractor 48 is trying to change the values of pixels in a non-existent 
applied to a second input terminal of multiplexer 50. M+ 1 row. To block changes to die gray-scale values of 
The output signal from comparator 46 is used to select the later occurring pixels where there are no later oc- 
bctwcen the two input signals to multiplexer 50, In 10 earring pixels in the row (x>N) and where there is no 
comparator 46, z^^ is compared to the threshold value, additional row (y>M), the signals from delays 106 
T. If zjcj^is less than T. then z,j,is selected as the output through 110 have to be switched to register 32 by corn- 
signal from multiplexer 50, and if Zx^ is greater than or parators 83 and 85. From I/O interface 22 the values of 
equal to T, then the output signal from subtractor 48 is N and x, and M and y are applied to comparators 83 and 
the output signal from multiplexer 50. The signal from 13 85. respectively. The output signal from each of corn- 
multiplexer 50 is Ex^ as discussed above. The output parators 83 and 85 are then inverted by mverters 87 and 
signal from comparator 46. z'xf, also represents the 88, respectively. The output signal from delay 96 must 
gray-scale value of the pixel at x.y that is to be printed, be blocked if either of the output signals from compara- 
displayed or stored by the present invention. Since the tors 83 and 85 is high, and the output signals from delays 
output signal from comparator 46 is a one-bit binary 20 98 and ^ must be blocked if the output sipal from 
signal, its value is cither 0 (pixel off) or I (pixel on). In comparators 83 and 85 is individuaUy high, rcspec- 
printing, a "pixel ofT' results in no color being applied tively. Switches 95 through 99 under the control of 
to the medium at that location, and thus a white pixel comparators 83 and 85 perform the necessary blocking 
remains (or whatever the color of the material on which when the corresponding condition is met by opcmng 
the printing is being performed), whereas on a monitor 23 the switch. At all other times the switches are closed, 
a "pixel ofr* results in a black pixel at that location. For Since switch 97 must be controlled by both compara- 
a "pixel on" signal the result is just the opposite. tors 83 and 85. the output signals from mverters 87 and 

The second path of this embodiment begins with the 89 are appUed to the input terminals of a two-input OR 

value of Zxj, as the address pointer to range function gate 91 to generate the combined control signal for 

look-up table 52 where 256 values of Wmin and Wmoxarc 30 switch 97. 

stored, one value of each for each potential value ofzx^. Second Embodiment 
This could also be implemented by a combination of r • • r 
subtractor. divider and multiplier to individually calcu- In the second embodiment of the present mvcnton of 
late the values for w„i>. and w««. The WmiV, is next an error diffusion process, codmg errors are randondy 
subtracted from w„„by subtractor 58 with the result of 35 diffused from each given pixel location to four immedi- 
that calculation being applied to one input terminals of atcly-adjaccnt ncighbonng pixel locations that are 
each of multipliers 68 and 70. The other input terminals scanned after the given pixel location. More particu- 
of each of multipliers 68 and 70 receive a random frac- larly. for each pixel location from which a codmg error 
tional value from look-up table 62. The output signals is diffused, a random choice is made of Uie order m 
from multipliers 68 and 70 represent the values n and r3, 40 which errors are to be diffused to the immediately-adja- 
respcctivcly. The value of ri is then siibtracted from 1.0 cent neighbors of the pixel locaUon under consider- 
by subtractor 72, and the value of rj is then subtracted ation; that is, the order of directions withm the direction 
from the resultant value from subtractor 72 by sub- set are randomly selected. For purposes of the foUow- 
tractor 74 to form the value of n. ing discussion, the directions from pixel P to the four 
The two initial paths of this embodiment now come 45 immediately-adjacent pixel neighbors that are scanned 
together at multipUers 76 through 80 where the value after pixel P are designated as directions di, d> dj, and 
Ex.v is multiplied by each of the values of ri, n and rs, d/ to emphasize that the directions arc randomly as- 
nKpectively; to form the values em. erri and crry The signed. To provide random assignments withm the di- 
ncxt step is to error correct the three later presented rection set, it is usually convcmcnt to use a sequencmg 
pixel values to which vectors di.dzanddj point in FIG. 50 look-up table containing approximately several thou- 
2. This is done by adders «« through 92 where the val- sand random numbers. 

ues of cm through erra arc added to z<x+ij.> 2<x-Hj;+l> Further in the second embodiment of an error diffu- 

and X(xj+\)* respectively, to form the new values of sioh process according to the present mvention, ancr 

those pUels as influenced by the coding error associated direction assignments are randomly determined withm 
with pixel Px Keep in mind that each pixel is influ- 53 a direction set, a first direction d/is assigned a weighting 

enccd by the error at up to three otiicr pixels in this factor Wrf/ where the value of the weighting factor is 

embodiment of the present invention. Finally, these randomly chosen from a uniform distribution of weights 

new values of those pixels are m^ed to delays 94 between 0.0 and 1.0. Next, a second randomly chosen 

through 98 with the output signals from the delays direction dj is assigned a weighting factor w^y where the 
being applied to register 32 to substitute the new values 60 value of the weighting factor is randomly chosen froin 

for the previous values contained in register 32 which a uniform distribution of weights between 0.0 and 

are shown on the associated output lines as z'. The (LO-wrfd- Next, the third randomly selected direction 

labels on the output lines from delays 106 tiirough 110 d*is assigned a weighting factor wrf* where the value of 

are z'V, i Di z' W+ W) which correspond to the weighting factor is randomly chosen from a uniform 
their locations after the occurrence of the next byte 65 distribution of weights between 0,0 and 

clock signal. (1.0-Wrf,-w<(/). FinaUy, the remaining direction d/ is 

The fmal piece of housekeeping is to repeat the above assigned a weighting factor wj/ such that 

for each pixel in the image in response to the byte clock wdi^ 1 - wdi- w^f- wok^ Thus, it may be noted that: 
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pixels located at (x+5y,y+5), (x+5y,y+6), 
1- (x+4,y+6), and (x+3,y+6). 

For multi-color halftone printing, a field such as that 
As mentioned above, pseudo-random weights vidi, jho^ „, piQ. 2 can be sequentially scanned with dif- 
wrf> Wrf*. and wji are used to propagate, or spread, 5 color filters or, alternatively, can be scanned 

coding errors to the four immediately-adjacent neigh- ^j^^ each picture element having multiple expo- 
bors of each scanned pixel location. For convenience of ,urcs with different filters. Then for each of the selected 
expression, the propagated error values are designated ^ ^^^^ pj^^j locations in the scanned 

herein as ei, ej, ej. and e4 to designate the error valuw jj^j^ „g ^^^^ ^ described 

which are distnbuted to pixd locations (x+l.y), ,0 above. Multi-color images without discernible artifacts 
(x+ \.y+ 1). (x,y-H). and (x- l,y + 1). /f P^^vely. ^Oitnbt printed by conventiomd digital halftone 
Also, for convemence It can be assumed that weighting orintinB Brochures 

factor w* is associated with error e,, that weighting P For dLr monochrome or multiK»lor halftone print- 
factor w,j IS associated with error ej, that weighnng " above-described process can be summariz^ by 
factor wdk is associated with error ea, and that ,< »c noovc-orav-now vou »imuuai«,c« 1,7 

weighting factor wji is associated with error e*. the following algorithm: 

According to the foregoing, the error value e, which . »> ^°L'^^'^ P«=' T't^rj^-^ P«f' Z"^' 
is to be propagated to tht pUel located at (x-h l.y) from fo". value E, of the ith puel iis follows: 

the pixel located at (x,y) is determined by multiplying .=fy,.i£ 
the weighting factor wjjby the encoding error for pixel jo 

Tit^f l^;ifi:/^i''&2i !i^t'!^n*^r^Tw'Jvk°''V where e,is the error value diffused to the jth pixd loca- 

gated to the pixel located at (x,y+l) IS e2=(wrf»X*^x^^ ^ £_ • i i ^ i. ux. • i 

Ld the error ea propagated to the pixel loc'ated^at fr^"* f P«^l location where &e jth puel 

(x+l,y+l) is e3=(wd;t)(Ex^). Finally, the error e4 ^^^Jf" ^ V^!^^fy 1^^'^!^^^,^?^ the predeter- 
propagated from pixel P to the pixel located at „ mined^set of pucls neighbonng the ith pixd^ 
(x-i,y+l) is eAH^diKB,.) CoUectively, the error " and whwewyis the wcightmg factor determmed for t^^ 
values ei, ez, es, and C4 can be understood to be pseudo- J* location according to the procedures descnbcd 
randomized error values. It should be noted that those above; 

propagated error values can be either positive or ncga- 2) continue the error diffusion process for the other 
^ve. P^^' neighbors of the ith pixel location; and 

As a result of a coding error being diffused, gray- 3) continue the error diffusion process for the (i+ l)th 
scale values of neighboring pixels are increased, or de- pixel location, and so forth. 

creased. For example, after error diffusion from pixel P FIG. 6 shows _a block diagram implementation of 
at location (x,y), the gray-scale value imputed to the computer 7 for the second embodiment enor diffusion 
location (x-hUy) would be the detected gray-scale techiuque of the present invention. This implementation 
value for that location plus the propagated error ei. is for a pixel pattern as shown in FIO. 2 with each of the 
Likewise, the gray-scale value imputed to the pixel registers initially set to zero, N pixels to a row, and M 
located at (x + liy + 1) would be the detected gray-scale pixel rows. The z gray-scale vdues for each pixel are 
value for that pixel location plus the propagated error sequentially applied to an N+1 byte>long register 32 
C2. The effects of the error diffusion process on the pixel from I/O interface 22 that contain the z values for pixels 
neighbors of pixel P in this example can be summarized p^^^ through P(xo'+i> ^th the z value for P(x+l^+i) 
as follows: being the next data to be received by register 32. From 

register 32 the z values of each of the five pixels needed 
piMi (jr-no')=pi«>u+i^)+'i for the operation of the second embodiment of the pres- 

43 est invention are available. 
. . . , There are two paths in this operation which come . 

puei „-p«ei<,^.,^+ 1)+^ together at the end for the final calculations. First, there 

is the path for calculating the error associated with the 
pixei'(;v+i)=pUd(jv+i)+e3 and Operational pixel, Also, Zjcjp is applied to compara- 

3Q tor 46, the positive terminal of subtractor 48, and multi- 
plexer 50. In subtractor 48 the value 255 is subtracted 
pixei'(jr.iji.i)«pixei(x.i^+i)+«4 ftom Zx^ and the resulting value from subtractor 48 is 

applied to a second input terminal of multiplexer 50. 
where the primes indicate that the gray-scale values of The output signal from comparator 46 is used to select 
the pixel locations are updated as a result of the above- 55 between the two input signals to multiplexer 50. In 
described error diffusion process. comparator 46, Zx^ is compared to the threshold value, 

It should be emphasized that the above-described T, If Z;^^ is less than T, then Zx^^ is selected as the output 
error diffusion process continues for each successively signal from multiplexer 50, and if Zxj^is greater than or 
scanned pad location* For example, as a result of the equal to T, then the output signal from subtractor 48 is 
error diffusion process, the gray-scale value for the ^ the output signal from multiplexer 50. The signal from 
pixel located at (x-l-4,y-|-5) will reflect the error propa- multiplexer 50 is Ejk^ as discussed above. The output 
gated from its neighboring pixels located at (x+3,y-f.4), agnal from comparator 46, z'",j^ also represents the 
(x+4,y 4-4), (x-f-5y.y-l-4), and (x-f 3,y-l-4). Also, the gray-scale value of the pixel at x,y that is to be printed, 
binary coding error propagated from the pixel location displayed or stored by the present invention. 
(x+4,y-H5) will reflect the errors propagated to that ^ The second path of this embodiment begins with the 
location from its set of neighboring pixels, and that presentation of three random numbers, each having a 
error, in turn, will be weighted by the above-described fractional value from look-up table 200. One of these 
pseudo-random process and then will be diffused to the random numbers is arbitrarily selected for wdghting 
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factor wdi. That value is then applied to subtractor 222 
where it is subtracted from 1.0 and the result applied to 
multiplier 208. A second of the random number from 
look-up table 200 is also applied to multiplier 208 to 
generate the second wdghting factor wjy. Weighting 5 
factor v/dj is then subtracted from the resultant value 
from subtractor 222 by subtractor 230. The resultant 
value from subtractor 230 is then multiplied by the third 
random number from look-up table 200 to generate 
weighting factor wrfjt. Ncxt» weighting factor wja is 10 
subtracted from the resultant value from subtractor 230 
by subtractor 238 to form weighting factor Wrf/. 

At this point the two paths merge with the individual 
multiplication of each of the weighting factors by the 
error signal Exj^in multipliers 242 through 248 with the 15 
resultant values being error values c/, cy, ejk and ei. The 
next step is to randomly collate these error values with 
each of the pixels to which vectors di through d4 are 
pointing in FIG. 2. To do that, four multiplexers 256 
through 262 have been included with each of c/ through 20 
e/ being applied as the four input signals to each of 
multiplexers 256 through 262. The selection from the 
input signals is performed imder the control of look-up 
table 252 which contains, in random order, all of the 
possible combinations of the four numbers 1 through 4. 25 
Thus, the errors values are provided by multipliers 256 
through 262 as ei through 64, respectively. The next 
step is to error correct the four later pixel values to 
which vecton dj, 6h da and d4 point in FIG. 2. This is 
done by adders 264 through 270 where the values of ei 30 
through C4 arc added to z(x+ 1 1 1> z(xjr+ 1) and 
Z(x-}^+ 1), respectively, to form the new values of those 
pixels as influenced by the coding error at Keep in 
mind that each pixel is influenced by the coding error at 
up to four other pixels in this embodiment of the present 33 
invention. Finally, these new values of those pixels are 
applied to delays 272 through 278 with the output sig- 
nals from the delays being applied to register 32 to 
substitute the new values for the previous values con- 
tained in register 32. 40 

The final piece of housekeeping is to repeat the above 
for each pixel in the image in response to the byte clock 
signal, and when x equals 0 or N or y equals M, to block 
the changes to the pixel values in register 32 from delay 
278 (x=0). delays 272 and 274 (x«N) and delays 45 
274-278 {y=M), respectively. The blocking discussed 
above is necessary to prevent the value of the pixels in 
the first column from effecting the pixels in the last 
column (x=0), the pixels in the last column from effect- 
ing the pixels in the first column (xs=N)» and from try- 50 
ing to change the values of pixels in a non-existent M + 1 
row. To block changes to the gray-scale values of the 
later occurring pixels where there are no earlier occur- 
ring pixels in the row (x< 1), where there arc no later 
occurring pixels in the row (x>N) and where there is 55 
no additional row (y>M), the signals from delays Z72 
through 278 have to be switched to register 32. 

The circuitry for performing this function is shown in 
FIG. 7. The values of x, y, N and M are obtained form 
I/O interface 22 in FIG, 4, and signals A through D are 60 
obtained from FIG. 6. Comparators 300 through 304 
determine whether x=0, x=N and y=M. respectively. 
The output signals from each of comparaton 302 and 
304 are then inverted by inverters 306 and 308, respec- 
tively. The signals from inverters 306 and 308 individu- 65 
ally control the switching of switches 320 and 318, 
respectively, to block B and D when the appropriate 
condition is met. Signals A and C need to be blocked 
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when either of two conditions are met; these switches 
314 and 316 are individually controlled by the output 
signal from AND gate 312 and 310, respectively. OR 
gate 316 allows control of switch 314 by cither compar- 
ator 300 or 304, and OR gate 312 allows control of 
switch 316 by either comparator 302 or 304. 

In conclusion, although the preferred embodiments 
of the present invention have been described, those 
s^led in the art will appreciate that additions, modifi- 
cations, substitutions and deletions not specifically de- 
scribed in the foregoing may be made to the preferred 
emt>odiment of the process. For example, diffusion di- 
rections can be chosen other than the above-discussed 
set of four directions. Moreover, it is not necessary for 
the present invention that errors be diffused to only 
immediately-adjacent pixel neighbors. 

Accordingly, the above-described embodiments 
should be regarded as illustrative rather than restrictive, 
and it should be appreciated that workers skilled in the 
art may make variations in those embodiments without 
departing from the spirit or scope of the present inven- 
tion as defined by the following claims. 
What is claimed is: 

1. An image enhancement process with error diffu- 
sion, comprising the steps of: 
detecting gray-scale values at preselected pixel loca- 
tions in an image; 
for each detected gray-scale value at the preselected 
pixel locations, determuiing upper and lower limits 
of a range fiinction which depends on the detected 
gray-scale values, and selecting first and second 
values which are randomly located between the 
upper and lower limits of ibt range function, and 
assigning a third error value for the preselected 
location equal to its detected gray-scaJe value mul- 
tiplied by 1 minus the sum of the first of the se- 
lected randomly located values and the second of 
the selected randomly located values; 
if the detected gray-scale value at a preselected pixel 
location is less than a predetermined threshold 
value, encoding the pixel location as a binary zero, 
assigning a first error value for the preselected 
pixel location equal to its detected gray-scale value 
multiplied by the first of the selected randomly 
located values and assigning a second error value 
for the preselected pixel location equal to its de- 
tected gray-scale value multiplied by the second of 
the selected randomly located values; 
if the detected gray-scale value at a preselected pixel 
location is greater than the threshold value, encod- 
ing the pixel location as a binary one, assigning a 
first error value for the preselected pixel location 
which equals the first of the selected randomly 
located values multiplied by ia quantity equal to the 
detected gray-scale value less 255, and assigning a 
second error value for the preselected pixel loca- 
tion which equals the second of the selected ran- 
domly located values multiplied by a quantity 
equal to the detected gray-scale value less 255 and 
assigning a third error value for the preselected 
pixel locations eqiial to its detected gray-scale 
value less 255 multiplied by 1 minus the sum of the 
first of the selected randomly located values and 
the second of the selected randomly located values; 
then, for each preselected pixel location, distributing 
the first and second assigned error values to at least 
three adjacent preselected pixel locations, one of 
which lies on the same line as the preselected pixel 
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location and the other two of which are on an 
adjacent line; and 
retaining images for display on a selected medium 
based upon values of the preselected pixel locations 
that have been encoded after distribution of the 
assigned error values. 
2. A digital halftoning process with random error 
diffusion for reducing visually apparent artifacts in im- 
ages that are produced by digital halftone image en- 
hancement» comprising the steps of: ^ 
detecting intensity values at selected pixel locations 
in an image; 

encoding selected ones of the pixel locations as one 
binary value when the detected intensity value at a 
location is less than a predetermined threshold IS 
value, and calculating error values for such pixel 
locations that equal the detected intensity values; 

encoding selected other ones of the pixel locations as 
a second binary value when the detected intensity 
values at those locations are greater than the prede- 20 
termined threshold value, and calculating error 
values for those locations that equal the detected 
intensity values less the upper limit of the intensity 
scale; 

for selected pixel neighbors of the ith pixel in an 25 
image* diffusing the calculated error value, E/, of 
the ith pixel location to its selected neighbors as 
follows: 



30 



where e/ts the error value diffused to the jth pixel 
neighbor from the ith pixel location, where the jth 
pixel neighbor is randomly determined from a pre- 
determined set of pixels neighboring the ith pixel, 35 
and where v/j is a weighting factor detennined for 
the jth pixel location; 

continuing the error diffusion process for the other 
pixel neighbors of the ith pixel location; 

continuing the error diffusion process for the (i -h l)th 40 
pixel location and so forth, such that coding errors 
are diffused from substantially all of the pixel loca- 
tions in a scanned image; and 

retaining images for display on a selected medium by 
digital halftone printing based upon intensity val- 45 
ues of pixel locations that have been modified by 
error diffusion. 

3. A digital halftoning process with random error 
diffusion for reducing visually apparent artifacts in im- 
ages that are produced by digital halftone image en- so 
hancement, comprising the steps of: 

detecting intensity values at selected pixel locations 
in an image; 

for a fint selected pixel location, encoding the pixel 
location as one binary value when the detected ss 
intensity value at the location is less than a prede- 
termined threshold value, and assigning an error 
value to the pixel location equal to the detected 
intensity value; 

when the detected intensity value at the first selected 60 
pixel location is greater than the threshold value, 
encoding the pixel location as a second binary 
value and calculating an error value for that loca- 
tion that equals the detected intensity value less the 
upper limit of the intensity scale; 63 

randomly choosing a first neighbor to the first pixel 
location from a predetermined set of neighboring 
pixels; 



for the first randomly chosen neighbor, assigning a 
weighting factor which is randomly selected; 

randomly choosing other neighbors to the first pixel 
location from the first predetermined set, and as- 
signing weighting factors to those locations that 
are randomly selected; 

for the selected pixel neighbors, diffusing error values 
as follows: 

where eyis the error value diffused to the jth pixel 
neighbor from the ith pixel where the jth pixel 
neighbor is randomly determined from a predeter- 
mined set of pixels neighboring the ith pixel, w; is 
the weight determined for the jth pixel location, 
and Ei is the error value assigned to the ith pixel 
location; 

assigning revised gray-scale values to the selected 
neighbor pixels such that the revised gray-scale 
values equal the detected gray-scale values plus the 
diffused gray-scale values; and 

retaining images for display on a selected medium by 
digital halftone image enhancement based upon 
intensity values of pixel locations that have been 
modified by error diffusion. 

4. Apparatus for image enhancement with error diffu- 
sion comprising: 

means for detecting gray-scale values at preselected 
pixel locations in an image; 

means for determining upper and lower limits of a 
range function which depends on the detected 
gray-scale values for each detected gray-scale 
value at the preselected pixel locations; 

means for selecting first and second values which are 
randomly located between the upper and lower 
limits of the range function for each detected gray- 
scale value at the preselected pixel locations, and 
assigning a third error value for the preselected 
pixel location equal to its detected gray-scale value 
multiplied by 1 minus the sum of the first of the 
selected randomly located values and the second of 
the selected randomly located values; 

means for determining if the detected gray-scale 
value at a preselected pixel location is less than a 
predetermined threshold value; 

means for encoding the pixel location as a binary rero 
if the detected gray-scale value at a preselected 
pixel location is less than a predetermined thresh- 
old value; 

means for assigning a first error value for the prese- 
lected pixel location equal to its detected gray- 
scale value multiplied by the first of the selected 
randomly located values if the detected gray-scale 
value at a preselected pixel location is less than a 
predetermined threshold value; 

means for assigning a second error vidue for the pre- 
selected pixel location equal to its detected gray- 
scale value multiplied by the second of the selected 
randomly located values if the detected gray-scale 
value at a preselected pixel location is less than a 
predetermined threshold value; 

means for determining if the detected gray-scale 
value at a preselected pixel location is greater than 
the threshold value; 

means for encoding the pixel location as a binary one 
if the detected gray-scale value at a preselected 
pixel location is greater than the threshold value; 
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means for assigning a first error value for the prese- 
lected pixel location which equals the first of the 
selected randomly located values multiplied by a 
quantity equal to the detected gray-scale value less 
255 if the detected gray-scale value at a preselected 
pixel location is greater than the threshold value; 

means for assigning a second enor value for the pre- 
selected pixel location which equals the second of 
the selected randomly located values multiplied by 
a quantity equal to the detected gray-scale value 
less 255, and assigning a third error value for the 
preselected pixel location equal to its detected 
gray-scale value less 255 multiplied by 1 minus the 
sum of the first of the selected randomly located 
values and the second of the selected randomly 
located values; 

if the detected gray-scale value at a preselected pixel 
location is greater than the threshold value; 

means for distributing the first and second assigned 20 
error values to at least three adjacent preselected 
pixel locations, one of which lies on the same line as 
the preselected pixel location and the other two of 
which are on an adjacent line for each preselected 
pixel location; and 

means for retaining images for display on a selected 
medium based upon values of the preselected pixel 
locations that have been encoded after distribution 
of the assigned error values. 

5, An apparatus for digital halftoning with random 
error diffusion to reduce visually apparent artifacts in 
images that are produced by digital halftone image 
enhancement comprising: 

means for detecting intensity values at selected pixel 35 
locations in an image; ' 

means for encoding selected ones of the pixel loca- 
tions as one binary value when the detected inten- 
sity value at a location is less than a predetermined 
threshold value; 

means for calculating error values for such pixel loca- 
tions that equal the detected intensity values; 

means for encoding selected other ones of the pixel 
locations as a second binary value when the de- 
tected intensity values at those locations are 
greater than the predetermined threshold value; 

means for calculating error values for those locations 
that equal the detected intensity values less the 
upper limit of the intensity scale; 

means for diffusing the calculated error value, Eh of 
the ith pixel location to its selected neighbors for 
selected pixel neighbors of the ith pixel in an image, 
as follows: 
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where eyis the error value diffused to the jth pixel 
neighbor from the ith pixel location, where the jth 
pixel neighbor is randomly determined from a pre- ^ 
determined set of pixels neighboring the ith pixel, 
and where wyis a weighting factor determined for 
the jth pixel location; 



means for continuing the error diffusion process for 
the other pixel neighbors of the ith pixel location; 
means for continuing the error diffusion process for 
the (i-l-l)th pixel location and so forth, such that 
coding errors are diffused from substantially all of 
the pixel locations in a scanned image; and 
means for retaining images for display on a selected 
medium by digital halftone image enhancement 
based upon intensity values of pixel locations that 
have been modified by error diffusion. 
6. An apparatus for digital halftoning image enhance- 
ment with random error diffusion for reducing visually 
apparent artifacts in unages that are produced by digital 
halftone image enhancement comprising: 
means for detecting intensity values at selected pixel 

locations in an image; 
means for encoding the pixel location as one binary 
value when the detected intensity value at the loca- 
tion is less than a predetermined threshold value 
for a first selected pixel location; 
means for assigning an error value to the pixel loca- 
tion equal to the detected intensity value for a first 
selected pixel location; 
means for determining when the detected intensity 
value at the first selected pixel location is greater 
than the threshold value; 
means for encoding the pixel location as a second 
binary value and calculating an error value for that 
location that equals the detected intensity value 
less the upper limit of the intensity scale when the 
detected intensity value at the first selected pixel 
location is greater than the threshold value; 
means for randomly choosing a first neighbor to the 
first pixel location from a predetermined set of 
neighboring pixels; 
means for assigning a weighting factor which is ran- 
domly selected for the first randomly chosen 
neighbor; 

means for randomly choosing other neighbors to the 
first pixel location from the first predetermined set; 
means for assigning weighting factors to those loca- 
tions that are randomly chosen; 
means for diffusing error values for the selected pixel 
neighbors, as follows: 

where cjis the error value diffused to the jth pixel 
neighbor from the ith pixel where the jth pixel 
neighbor is randomly determined from a predeteir* 
mined set of pixels neighboring the ith pixel, w/is 
the weight determined for the jth pixel location, 
and Ei is the error value assigned to the ith pixel 
location; 

means for assigning revised gray-scale values to the 
selected neighbor pixels such that the revised gray- 
scale values equal the detected gray-scale values 
plus the diffused gray-scale values; and 
means for retaining images for display on a selected 
medium by digital halftone image enhancement 
based upon intensity values of pixel locations that 
have been modified by error diffusion. 
• * • « * 
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